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Abstract 
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^SJ , We study the effects of small chemical potentials associated with the three light quark flavors in QCD. We use 

a low-energy effective field theory that solely relies on the symmetries of the QCD partition function. We find 
three different phases: a normal phase, a pion superfluid phase and a kaon superfluid phase. The two superfluid 
phases are separated by a first order phase transition, whereas the normal phase and either of the superfluid 
phases are separated by a second order phase transition. We compute the quark-antiquark condensate, the pion 

CN ' 

• condensate and the kaon condensate in each phase, as well as the isospin density, the strangeness density, and the 

c 2 ^ , mass spectrum. 
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a ' 1 Introduction 



QCD at non-zero baryon density has recently been the subject of many studies. Our present 
knowledge of nucleon rich matter is limited to low densities and very high densities. In the first 
regime, we can rely on our good phenomenological understanding of nuclear interactions in the 
vacuum. In the second regime, the phenomenon of color superconductivity and its consequences 
are becoming better understood both qualitatively and quantitatively [Q, ||, ^, However we 
have only a poor grasp of a sample of nuclear matter that is neither too dilute nor too dense. 
Unfortunately this kind of matter is found in many interesting physical systems, from supernova 
explosions and neutron stars to heavy ion collisions. In other words, our knowledge of the QCD 
phase diagram at non-zero chemical potential is limited to extreme cases, and it is clear that we 
need a better understanding of its characteristics. 

This situation is very different from what we know about QCD at non-zero temperature 
and zero chemical potential, where numerous analytic and numerical studies have enabled us to 
achieve a qualitative and a quantitative understanding . Unfortunately, the current numerical 
techniques are useless in the case of QCD with three colors and fermions in the fundamental 



representation at non-zero baryon density. The complexity of the measure in the partition 
function is the origin of the problem. Special techniques have been developed to tackle this 
problem, but their reach is presently far from QCD @. There are however QCD-like theories 
at non-zero baryon density that do not have this complexity problem. Therefore they can be 
studied numerically with the usual techniques: QCD with two colors and fundamental fermions, 
and QCD with any number of colors and adjoint fermions. We now have a good understanding 
of these QCD-like theories from analytic and numerical studies ||^ ^, ^, |l^, |l^, QCD 
with three colors and two flavors of fundamental fermions at non-zero isospin density also has a 
real measure, and can therefore also be studied with standard numerical algorithms [p!^]. 

In this article we will study the QCD phase diagram for small chemical potentials associated 
to each quark flavor. At low enough energies, the single most important property of QCD is that 
chiral symmetry is spontaneously broken in the vacuum. It implies the existence of light bosonic 
degrees of freedom, and stringent constraints on the physical observables. This property has 
been extensively and successfully used to study QCD at low energy within an effective theory 
solely based on the symmetries of the QCD partition function, Chiral Perturbation Theory 
p5| , |l^] . Since the spontaneous breaking of chiral symmetry is still the most relevant property of 
QCD at low enough quark chemical potentials, Chiral Perturbation Theory should also enable 
us to explore the QCD phase diagram in this domain. This approach has been successfully used 
for QCD-like theories at non-zero baryon chemical potential and for QCD at non-zero isospin 
chemical potential ||^, |8|, |9|, [m . 

We will restrict ourselves to the three light quark flavors, for which this effective theory 
approach can be used. There is however an intrinsic limitation to this approach. Since only the 
octet of Goldstone bosons are kept as the relevant degrees of freedom, high chemical potentials 
where other excitations become relevant cannot be reached. A rapid glance at the hadron 
spectrum, and in particular the proton, the rho and the omega, shows that we will be able to 
explore the domain restricted by \fj,u\ < 300MeV, < SOOMeV, and \fis\ < 550MeV. 

Our study extends the analysis of QCD at low isospin chemical potential that included the up 
and down quark flavors only [p!^ ]. In this analysis, it was found that the pions form a condensate 
at high enough isospin chemical potential. Similarly, we will find that for high enough strange 
chemical potential, it is energetically favorable for the kaons to form a condensate. A non-trivial 
phase diagram emerges; there are three distinct phases: a normal phase, a pion condensation 
phase and a kaon condensation phase. The superfluid phases are separated from each other by 
a first order phase transition, whereas either of them is separated from the normal phase by a 
second order phase transition. 

The pion and kaon condensates we will find in our study are very different from the pion 
and kaon condensates that emerge in a nucleon rich environment |17]. In this latter case, the 
condensation is driven by the attractive meson-nucleon interactions. In our case the baryon 
density is zero, and the pion and kaon condensates emerge out of the vacuum as soon as it 
becomes energetically favorable for the system to form such condensates in the ground state, 
just because of the charge carried by these modes. 

We proceed as follow. First we shortly derive the effective Lagrangian in Section 2, and 
then determine the ground state of the effective theory in Section 3. Section 4 is devoted to 
the computation of the various condensates relevant to the three different phases as well as the 
isospin and strangeness densities. The mass spectrum is computed in Section 5. Concluding 
remarks and discussion are presented in the Section 6. 
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2 Low-energy Effective Lagrangian 



In order to construct the low-energy effective Lagrangian, we need to study the symmetries of 
QCD. We work in Minkowski space. The fermionic part of the QCD Lagrangian at non-zero 
quark chemical potentials is given by 



where 



c 



QCD 




M = diag (m„,mrf,ms 



(1) 



(2) 



and B,y 



-B,0), B = diag{fiu,fJ.d,fJ's 



, A 11 11 ^ / ^ 

diag(-/iij + -/i/, -fiB - ^f^l, ^fJ-B - fJ-s)- (3) 



One can either use chemical potentials associated with the quark flavors, or chemical potentials 
associated with baryon number, ^b, isospin, fij, and strangeness, ^s- In this study, we will take 
rriu = rud = m. When all quark masses and chemical potentials vanish, the QCD Lagrangian 
is invariant under SUii^) x SUji{3) x [7^(1) • At zero chemical potential, this symmetry is 
spontaneously broken to SUv{3) x C/y(l). There are therefore eight Goldstone bosons in the 
spectrum: vr'', vr"*", vr", r/^, K~^, K~ , K^, and K^. They are the relevant degrees of freedom at 
low energy. 

The effective Lagrangian of QCD at low energies and finite quark chemical potentials in 
Minkowski space is given by 



F' 



cS 



^ TrV^^SV^St + ^GTtM{J: + 5]t), 



(4) 



where S G SU{3), and V,yS = 



Lagrangian at lowest order [15, 16 1 
spontaneous chiral symmetry breaking 



i[B^,Ti]. This is the usual Chiral Perturbation Theory 
It contains only the octet of Goldstone bosons due to 
The chemical potential has been introduced in the 

Under a 



same way as a regular vector source using the usual flavor gauge symmetry |16, | 
chiral rotation, the quark fields transform as {ipi, i'R) — > {VitpL, VRi^n), with Vl,r G SU{3)l,r, 
whereas S VrSV^. For zero chemical potentials S is invariant under vector transformations 
and contains the octet of Goldstone bosons due to spontaneous chiral symmetry breaking: 



S = U^U, with U = exp{i<^/V2F), 



and 



+ 



V 







\ 



V6 

K- 



vr ' 

^0 



K+ 

2r?" 
V6 / 



(5) 



(6) 



with S = diag(l, 1, 1) for ^b = /U/ = /U5 = 0. The effective Lagrangian is invariant under local 
chiral transformations. At zero chemical potentials, the masses of the Goldstone bosons are 
given by the Gell-Mann-Oakes-Renner relation: 



2 



2Gm/F^ 
G{m + ms)/F'^ 
2G{m + 2ms)/3F'^ 



{Am? 



K 



ml)/3. 



(7) 
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The first remarkable property of the effective Lagrangian @ is that in the {fiB, fJ-i, fJ-s)- 
basis, hb completely drops out of Ceg. This is just due to the fact that the effective Lagrangian 
contains only the octet of Goldstone bosons due to spontaneous chiral symmetry breaking. 
The use of the effective Lagrangian is therefore hmited to l/isl < 940MeV, < 770MeV, 
and lusl < 550MeV. These estimated bounds come from the masses of the proton, the rho 
and the omega, respectively. This means that the effective Lagrangian is only valid in the 
domain \nu\ < SOOMeV, < SOOMeV, and \fis\ < 550MeV. Therefore in this work, the 
baryon chemical potential will never appear explicitly, and we will work with and /x^. Recall 
however that all our results are valid for < 940MeV, and that Chiral Perturbation Theory, 
in its domain of validity, allows us to determine the phase diagram in the three quark chemical 
potentials. For definiteness, we will restrict ourselves to ni,fJ-s > 0, the other quadrants in the 
(/X/, /X5)-plane are easily derived from our results. 



3 Ground state 

The ground state of the effective theory is determined by the maximum of the static part of the 
effective Lagrangian. Therefore S, the maximum of 

Cst.t = -^Tr[B,S][5,St] + iGTr7W(S + St), (8) 

determines the ground state of QCD at non-zero quark chemical potentials. 
We will use the following Ansatz for the saddle point: 

/ 1 \ / cosa sina \ / 1 \ 

S = cos (3 — sin /? ~ sin a cos a cos /3 sin (3 \ , (9) 

^0 sin/3 cos/3 ) \ 1/ yO -sin/3 cos/3 J 

with a,/3 G (0, 7r/2). Notice that S|q,=o = diag(l,l,l) corresponds to the ground state at 
m = fj-s = 0- Furthermore in the massless case. 




S|(„=|,^=o) = I -1 I (10) 
is found to be the maximum when /i/ > 2iJ,s, and 




=1 10 1 (11) 

is found to be the maximum when /i/ < 2fis- By noticing that the condensate is basically given 
by 'i/'ijS'i/'L, it is easy to see that ( p^ corresponds to a pion condensation phase, and that (|ll]) 
corresponds to a kaon condensation phase. Our Ansatz corresponds to an axial rotation to go 
from the normal phase to either of the superfluid phases, and a vector rotation to go directly 
from one superfluid phase to the other. We will show that this Ansatz indeed corresponds to a 
local maximum of Cgtut- We will assume that it is the global maximum. 

After our Ansatz (P) is introduced into (§) , the static part of the effective Lagrangian becomes 
quite complicated 

-Cstat = ^^{(21f^j + l2^iI^s + 20^^l + {19fij + 20^^J^^s + l2^il)cosa)s■m\a/2) 
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-2(/xj - 2ns f sin^(a/2) cos(4/3) + 2{ni - 2/^5) (3/x/ + 2^x5) sin^ a} (12) 

+G|m cos a{l + cos^ /?) + cos^ /3 + (m + cos a) sin^ /jj . 

We first maximize Cgtat with respect to /3. Remember that we limit ourselves to the first quadrant 
in the (/x/, /i5)-plane: fii,fJ,s > 0. We find that the maximum, for a 7^ and /3 G (0, 7r/2), is 
either at /3 = or at (3 = tt/2. Therefore there are two different phases according to the value 
of p. We look for the maximum of £stat separately in each case. 



i) P = 0: The static part of the effective Lagrangian, its maximum, and the saddle point are 

F2 



given by 



with 
and 



2Gm cos a + Gmg + sin^ a, 

cos a = 1 ioT fj,i < m^r 

coscK = (m7r//U/)^ for fj,i > m^r, 

/ cos a sin a \ 

— sin a cos a 
1 



V 



(13) 
(14) 

(15) 



ii) (3 = 7r/2: The static part of the effective Lagrangian, its maximum, and the saddle point 



are given by 



with 
and 



1 

Gm + G{m + rus) cos a + "^(2^^ A*^)^ 

{cos a = 1 for i/// + jjls < niK 

cos a = {rriK film + IJ'S))'^ for 2!^! I^S > rriK, 

I cos a sina \ 



v 



1 

sin a cos a j 



(16) 
(17) 

(18) 



Finally, in order to find the maximum of the static part of the effective Lagrangian in the 
whole quadrant > 0, we have to compare the value of £stat for a = 0, for /3 = and 

cos a = (m7r//U/)^, and for /3 = 7r/2 and cos a = {m.K/i^fJ'i + /^s))^- The orientation of the 
condensate for a given pair {fii^fis) is determined by the maximum of -Cstat at this point. We 
find three different phases: 



1. Normal phase: m < 171^^ and /X5 < uik — 

cos a = 1, any j3 £ (0, 7r/2) 



(19) 



2. Pion condensation 



phase: m > and /X5 < (-m^ + ^J {ml - /x|)2 + Amj^ftj^ /2/x/ 



cos a = 



1^1 



(20) 



3. Kaon condensation 



phase: us > mK-\lii and /X5 > (-m^ + ^ {ml - /xf )2 + 4mj^nfj /2iii 



cos a 



rriK 



, P = 7r/2. 



(21) 
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The curve /i^ = ^- 



-ml + 



fif)'^ + Amj^fi^ 1 /2^/ is a first order phase transition curve, 



because the pion and kaon condensates compete on this curve. The phase diagram in the 
first quadrant and its obvious extension to the whole (///, //5)-plane are given in Fig. 1 for 
m,r = 140MeV, and rrix = SOOMeV. They are very similar to the phase diagram of QCD-like 
theories at finite baryon and isospin densities Q. 
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Figure 1: Phase diagram in the (/x/, /i5)-plane. The solid lines are second order phase transitions, 
whereas the dotted curves are first order phase transitions. The intersection of the three phase transition 
curves on the left figure is at (/i/, ^s) = {TnT^^mK — ra-„/2). 



4 Condensates and charge densities 



We will show in the next section that S given by (PJl9|,20,21) is indeed a local maximum of the 
static part of the effective Lagrangian. We will assume that it is the global maximum. But first 
we investigate the properties of the ground state described by (|g|,19|, 20| , ^T|) in the three phases 
by looking at different observables. The quark-antiquark condensates, {qq) for q = u,d,s, the 
pion condensate, {d'j^u + h.c), the kaon condensate, {s^^u + h.c.) are obtained by inserting 
the appropriate sources into the effective Lagrangian [Q, or more directly by noticing that the 
vacuum expectation value of an operator bilinear in the quark fields is basically given by (t/^piJ^ipL) 
( p^JTl| ). The isospin density, n/ = dCcs/dm, and the strangeness density, ns = dC^s/dfis are 
easily obtained from the static part of the effective Lagrangian. 



1. Normal phase, fij < and fis < — hfJ-i- 



{uu) = {dd) = {ss) = G= {qq)o 
{d'y^u + h.c.) = (s75ti + h.c.) = 
nj = ns = 



(22) 
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2. Pion condensation 



phase, m > rriTr and us < (-^"'^ + \J ("ml - Z^/)^ + 4m|-/Lt|) /2/X7: 



(ss) = G = {qq)o 



{dj^u + h.c.) = Gsina = (gg)o^l - 'mi/fij 
{s^5U + h.c.) = 
ni = F'^iii{l - ml/iJij) 
, ns = 



(23) 



3. Kaon condensation phase, /j^s > rriK—IJ'i/'^ and /is > (— + \J (fn^ — A*/)^ + ^fn\:l^^^ 



{uu) = {ss) = Gcosa= {qq)Qm\/{\ii'\ + iisf 
{dd) = G= {qq)o 
{dj^u + h.c.) = 



(s75U + h.c.) = Gsina = ((?g)oy 1 - m%/{2fJ'i + l^s) 
ni = + ^5) (1 - m%/{y.j + f^sY) 

ns = F^{^fii + us) (1 - m%/{kl^i + l^sY 



(24) 
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Figure 2: The figure on top is the phase diagram in the first quadrant of the (/x/, /i5)-plane. The hght 
solid horizontal line corresponds to us = 200MeV, the value we use throughout the other figures depicted 
underneath. In these last six figures, the critical value of /ij at the second order phase transition is 
depicted by a solid vertical line. The different condensates are shown in the next five figures. The isospin 
density (solid curve) and the strangeness density (dashed curve) are shown together in the last figure. 
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The isospin and strangeness densities in the superfluid phases just above the second order phase 
transition curves can be reproduced within a semi-classical analysis of a dilute Bose gas in the 
Bose-Einstein condensation phase ||8|. In the kaon superfluid phase, the strangeness density is 
naturally found to be twice bigger than the isospin density, since a kaon carries twice more 
strangeness than isospin. 

The different condensates and charge densities as a function of and for a constant us 
are depicted in Fig. 2, for fis = 200MeV, and in Fig. 3, for fis = 460MeV. In the first case 
we encounter only two phases as fij changes, whereas in the second case we encounter the 
three different phases found above as /i/ changes. As expected, the condensates and charge 
densities are continuous across the second order phase transition and discontinuous across the 
first order one. The different phases and the nature of the different phase transitions are clearly 
distinguished by the observables we study. 
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fj,j (MeV) /*j (MeV) (MeV) 



Figure 3: The figure on top is the phase diagram in the first quadrant of the (/i/, /j,5)-plane. The 
light solid horizontal line corresponds to jjs — 460MeV, the value we use throughout the other figures 
depicted underneath. In these last six figures, the critical value of /x/ at the second order phase transition 
is depicted by a solid vertical line, and the critical value of fii at the first order phase transition is 
depicted by a dotted vertical line. The different condensates are shown in the next five figures. The 
isospin density (solid curve) and the strangeness density (dashed curve) are shown together in the last 
figure. The difference between the second and first order phase transitions are very clearly seen in these 
observables. 



5 Curvatures and mass spectrum 

In this section we will prove that S given by (^Jl9, 2^,21) is indeed a local maximum of the static 
part of the effective Lagrangian. We compute the masses of the different excitations in the three 
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phases. In order to obtain these masses, we have to expand the whole effective Lagrangian up 
to second order in the Goldstone fields. The general form of the effective Lagrangian reads 

C^s = -^Tr[5,S][S,i;t] + iGTrA^(S + St) 

+— Tr9^S5^St + i— Tr5[St,5oq. (25) 

The fourth term in (^) may involve a term quadratic in the Goldstone fields (^). In the 
superfluid phases, the original Goldstone fields are mixed Since the usual Goldstone manifold 
changes as the normal phase is left, the original octet of Goldstone bosons is not an appropriate 
set of coordinates for the new Goldstone manifold. We will denote these new appropriate fields by 
a "tilde" . We choose them so that the kinetic term in the Lagrangian is canonical. Furthermore, 
the linear derivative term in (^) mixes these new fields. Therefore the mass spectrum is not 
directly given by the curvatures of the static part of the Lagrangian at the maximum. When 
n-fields are mixed, the inverse propagator that appears in the effective Lagrangian is a n x n 
matrix, and the masses are given by the zeros of its determinant in the pg-plane. A secular 
equation has therefore to be solved. It is a n-th order equation in pg when n fields are mixed 
by the linear derivative term in (p5|). These equations can be explicitly solved, but the explicit 
solutions are cumbersome and rather useless expressions. Therefore we will only give which 
fields are mixed by the linear derivative term in (^5|), and solve the secular equation explicitly 
in the simple cases only. The labels of the different mixed states will be given according to the 
usual Goldstone mode they correspond to at the second order phase transition with the normal 
phase. We study each phase separately, and summarize our results for the mass spectrum in a 
few figures. 



5.1 Normal phase 

For < rn-TT and fis < ""^x ~ ^A*/? with the Goldstone fields defined in (|6|), we find that there is 
no mixing induced by the linear derivative term in the quadratic part of the effective Lagrangian. 
The various masses read 

m^o = m^o = \J {Am?^ - ml)/3 

771^^+ = — m 

m^- = m^ + m (26) 

mK+ = rriK - hfJ-i - l^s 

rriK- = rriK + + l^s 

rriKO = rriK + h^^i - IJ-s 

rrij^o = niK - 2/"/ + /^5- 

The effect of the chemical potentials is just to shift the mass of a state by its charges times the 
chemical potentials: E ^ E — Ifij — Sfis, where / and 5 are the isospin and strangeness of that 
state. 



5.2 Pion condensation phase 

For m > m,r and fis < + \J ("^^ ~ /^/)^ + /2/i/, we find that there is mixing 

between some of the fields defined in (P), namely (vr", 77", vr'^, vr"), and {K~^ , K~ , , K'^). We 
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find that the different masses are given by 



2 



+Hj{4.mj^ + njf - Aml^i'](2Am\ - 43;u|) - Aml^xj{Am\ + iJ,j) 


^(irnlfij - 2mlfij + 6mi + 7fij - [^SmUjij - ml) 
+^|(4m2^ + - 'iminj{24mj^ - 43/i|) - AmlfijiAmj^ + nj) 

~+ ~- ~n -° 

and the masses of , K , K , K are given by the zeros on the E'-plane of 
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where 



s\ 

S2 
S3 
S4 



\{-m\ - i/if + + /i/(i^7 + /is)cosa) 
i{^^s + ^M/ cos a) 
^/// sin a 

ii^s - cos a), 



and cos a = (m^/^u/)^, as given in (pC 



(27) 



(28) 



(29) 



5.3 Kaon condensation phase 



For Us > mx — fJ'i/'^ and fis > (^—ml + \J {ml — fij)'^ + Amj^njj /2/X7,we find that there is 
mixing between some of the fields defined in (^), namely (vr'^, rj^,K~^, K~), and (vr^, vr^ 
The masses of Tr^,fj^, K~^, K~ are given by the zeros in the i?-plane of 
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i(— cos a + {^muj + /is)^ cos 2a), 



(31) 



and cos a = {mK/{^fJ'i + Ms))^, as given in (p|). One of the modes, the K^, is massless. The 



- 



masses of fr^ , vf , K , K are given by the zeros on the i?-plane of 



det 



V 



Eul 








EU3 







-EU3 

Eu^ 

+ U4 y 



(32) 
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where 



ui = ^{-2ml + fij - 2^1 us + + fJ-s) cos a) 

U2 = |(3^/ - + 2(i^/ + ^5) cosa) 

U3 = 2 + Ms) sin a (33) 

U4 = I (— 4m|^ + 2fij — 4:iJ,jfj.s — cos a + cos a) 

U5 = ji-3fj,i + 2fis + 2{^ni + fis) cos a), 



and cos a = {rriK/i^fJ-i + Ms))^) given in (ETI 
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Figure 4: The masses of the eight pseudo-Goldstone bosons are shown here as a function of /x/ for 
Us = 200MeV. In these two figures, the critical value of /i/ at the second order phase transition is 
depicted by a solid vertical line. 



The spectrum is depicted as a function of /// for us = 200MeV in Fig. 4, and for ^s = 
460MeV in Fig. 5, the corresponding condensates and densities are given in Fig. 2, and in 
Fig. 3, respectively. The masses are continuous across the second order phase transition lines, 
as they should, and they are in general discontinuous across the first order phase transition 
curves, as they might. In each superfluid phase there is one single massless Goldstone mode. 
The SUv{3) X Uv{i) symmetry that leaves the quark-antiquark condensate invariant is further 
broken by the new condensates in the superfluid phases: Ui{l) symmetry is broken by the pion 
condensate, and Us{l) is broken by the kaon condensate. 



6 Conclusions 

In this article we have explored part of the phase diagram of QCD at non-zero quark chemical 
potentials and zero temperature. We have used a low-energy effective theory based on the 
symmetries of the QCD partition function, Chiral Perturbation Theory. The quark chemical 
potentials have been introduced into the effective Lagrangian as a usual vector external source 
via a flavor gauge symmetry. Because we only kept the octet of the Goldstone bosons as relevant 
degrees of freedom, which is only justified for small enough chemical potentials, our analysis is 
intrinsicly limited to the domain < 300MeV, |/Ud| < 300MeV, and < 550MeV. But this 
whole domain can be reached within Chiral Perturbation Theory. 

We have found a non-trivial phase diagram. There are three distinct phases: a normal phase, 
a pion superfluid phase, and a kaon superfluid phase. The two superfluid phases are separated 
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Figure 5: The masses of the eight pseudo-Goldstone bosons are shown here as a function of /i/ for 
/i5 — 460MeV. In these two figures, the critical value of /i/ at the second order phase transition is 
depicted by a solid vertical line, and the critical value of /i/ at the first order phase transition is depicted 
by a dotted vertical line. 



by a first order phase transition curve, whereas each of them is separated from the normal phase 
by a second order phase transition curve. The observables that we have analyzed show a typical 
behavior for such a phase diagram. The behavior of the different condensates at the second order 
phase transitions is very similar to what has been found in QCD-like theories The number 

densities just above the second order phase transition curves can be understood from a dilute 
gas approximation The excitation spectrum in the normal phase is trivial, and there is one 
massless mode in each superfluid phase that is due to further spontaneous symmetry breaking 
by the unusual condensates, similar to what as been observed in QCD-like theories 0, ^. 

Our study extends the analysis of QCD at non-zero isospin chemical potential done by Son 
and Stephanov [14|. We recover all their results on the pion condensate, the isospin number 
density and the mass spectrum for < mx — ^Tr/S = 430MeV, where kaon condensation can 
become energetically more favorable, depending on (cf. Fig.l, 2, 4). 

Finally, at low energy and within the effective theory approach we have used, it would be 
interesting to study matter that is not only made out of quarks and gluons, but also electrons 
and neutrinos. The inclusion of the electro-magnetic interaction and of the weak interaction into 
Chiral Perturbation Theory has already been achieved in the vacuum [^] , and its generalization 
to non-zero quark chemical potential seems to be rather straightforward. The effective theory 
technique we have used in this article is therefore very well suited to study more general problems 
of this kind. 
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